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ABSTRACT 

Cosmic rays in the energy range from about lO's GeV to several lOO's TeV are observed on Earth 
with an energy-dependent anisotropy of order 0.01-0.1%, and a consistent topology that appears to 
significantly change at higher energy. The nearest and most recent galactic cosmic ray sources might 
stochastically dominate the observation and possibly explain a change in orientation of the anisotropy 
as a function of energy. However, the diffusion approximation is not able to explain its non-dipolar 
structure and, in particular, the significant contribution of small angular scale features. Particle 
propagation within the mean free path in the local interstellar medium might have a major role in 
determining the properties of galactic cosmic rays, such as their arrival distribution. In particular, 
scattering on perturbations induced in the local interstellar magnetic field by the heliosphere wake, 
may cause a re-distribution of anisotropic cosmic rays below about 100 TeV toward the direction 
of the elongated heliotail and of the local interstellar magnetic field in the outer heliosphere. Such 
scattering processes are considered responsible of the observed TeV cosmic ray global anisotropy and 
its fine angular structure. 

Subject headings: Magnetic fields - Scattering - Turbulence - Sun:heliosphere - ISMxosmic rays 



1. INTRODUCTION 

Galactic cosmic rays possess a small but signif- 
icant anisotropy of order lO"** — 10"'^. Obser- 
vations in the northern hemisphere were reported 
from energies of t ens to several hundred s GeV with 
muon detectors ( iNagashi ma "eTall [19981 : IHall et al.l 
Il999t iMunakata et al. 201(|, in th e 1-10's TeV en- 
ergy range with Tibet AS7 array (Amenomori et al.' 
[2OO61 1201 lal). Super-Kamiok ande ((Guilhan et al. 2007 
Milag r o (lAbdo et all I2OO9I) and ARGO-YBJ dZhan 
Imm IShuwand \2(im . up to lOO's TeV by EAS- 
TOP (lAs-lietta et all [20M) . The IceCube Ob- 
servatory recently reported similar observations in 
the southern hemisphere in the 10-100 TeV energy 
range (jAbbasi et al.l 120101 |2012D. also using the IceTop 
surface array component (jSantander et al.l l2012j ). The 
first major feature of the cosmic ray anisotropy is that 
it is more complex than a dipole, and that its ampli- 
tude seemingly increases with energy up to the lO's TeV 
range, above which it decreases again (l Amenomori et ahl 
[200 5). The second interesting characteristics is a sig- 
nificant change in topology at energies in excess of ap- 
proximately 100 TeV (as observed by EAS-TOP and Ice- 
Cube). And the third major property of the observations 
is the discovery of significantly small angular features in 
the TeV energy range. 

The origin of such uneven global distribution in the 
arrival direction of galactic cosmic rays is still unknown. 
A dipole anisotropy would be expected from the mo- 
tion of the solar system with respect to the cosmic ray 
plasma rest frame. If cosmic rays do not co-rotate with 
the Galaxy, an observer on Earth would see an ap- 
parent excess of particle counts toward the direction 



of galactic rotation and a deficit in the opposite di- 
rection. Such a convective anisotropy is referred to as 
the Compton-G etting effect (jCompton fc Getting! 119351 : 
iGleeson fc Axfo rd 1968). The predicted amplitude is ex- 
pected to be of order 10"'^ with the apparent relative 
excess in right ascension at about 21 hr. Apart from the 
fact that the observed anisotropy is not a simple dipole, 
neither the expected amplitude nor the phase seem to be 
consistent with the observations. On the other hand, the 
reference frame of the galactic cosmic rays is not known, 
therefore it is reasonable to assume that the Compton- 
Getting effect could be one of the contributions to the 
cosmic ray anisotropy. 

It has been speculated that the observed anisotropy 
might be a natural consequence of the stochastic na- 
ture of cosmic ray galactic sources, in particular nearby 
and recent Supernova Remnants (SNR). Discreteness 
of such sources, along with cosmic ray propagation 
through a turbulent interstellar medium, might lead 
to significant fluctuations of their intensity in space 
and, the refore, to a faint broad aii is otropy (see for 
instance [Erlvkin fc Wolfendale (20061); IBlasi fc Amatol 
l|2012f ): |Pohl fc Eic hler (2012)). In absence of signatures 
of nearby sources of cosmic rays an anisotropy could be 
generated by convection from persistent cumulativ e mag - 
netized flow fleld of old SNRs (jBiermann et al.l 120121) . 
Or it can result from the combined effects of the regu- 
lar and the turbulent (fluctuating) magnetic fields in our 
vicinity ([Battaner et al.ll2009t ). However, the solar sys- 
tem is surrounded by a h ighly heterogeneou s Local Inter- 
steUar Medium (LISM) ( Frisch et aLll20lH) . that can af- 
fect the arrival direction distribution of the cosmic rays. 
The common velocity of the nearby interstellar clouds 
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could indicate that they are part of an evolved sub-shell 
of the superbubble associated to the Loop I, expanding 
from the Scorpion-Centaurus Association. Various in- 
direct determinations of the Local Interstellar Magnetic 
Field (LIMF) suggest that its direction is coherent over 
scales of about 100 pc and roughly parallel to th e local 
surfac e of Loop I shell dFris ch 2010, .2011; Fris ch et alJ 
l2012al| bl). The observations, therefore, could be associ- 
ated to cosmic ray diffusive streaming between collid- 
ing n earby interstellar cloud s of the expanding Loop I 
shell (jSchwadron et al] 120121) . Since the solar system is 
located almost at the edge of the so-called Local Inter- 
stellar Cloud (Lie), a partially ioniz ed cloudlet within 
the Local Bu bble, it was proposed by lAmenomori et al.l 
((2007ll2011bD that a non-dipolar anisotropy could be gen- 
erated by the diffusion of cosmic rays through the LIMF 
connecting the solar system to the interstellar medium 
outside the LIC. On the other hand, the observation of 
a topological change of the anisotropy pattern at an en- 
ergy in excess of about 100 TeV, where a relative deficit 
is observed in the region of the sky dominated by a broad 
excess at lower energies, is an indication of a phenomeno- 
logical transition in the cause of high energy cosmic rays 
anisotropy. A proton in a 3 fiG magnetic field has a max- 
imum gyro-radius Rg ~ 80- ExeV AtQ. Even though the 
mean free path of cosm ic rays in magnet ized plasma can 
be significantly larger ()Yan fc Lazarian| [2b08V it is rea- 
sonable to believe that cosmic rays with energy below 
about 100 TeV must be affected by the heliosphere, con- 
sider ing its extension from hu nd reds to sev eral thousands 
AU ( iPogorelov et al.l l2009airbl: Izmodenov fc Al exashov. 
12003': Izm odenov fc Kallenbachll2006f r 



The anisotropy of TeV cosmic rays shows evidence of 
statistically significant small angular features. By effec- 
tively averaging the cosmic ray counts over an angular 
range of a few tens degrees, and subtracting the aver- 
aged map from the observation, the regions where cos- 
mic ray intensity changes significantly within such an- 
gular scale are identified as localized fractional excess 
or deficit regions. Such technique led to the discovery 
of two highly significant localized fraction al excess re- 
gions in the northern hemisphere by M ilagro (lAbdo et al 
2008f). also observed by Tibet AS7 (lAm enomo ri et al 



2006) and ARGO-YBJ (iVernetto et al.l200 9: lupp a et al 
20 Hi) . Similar observations were recently reported 



in the souther n hemisphere by t he IceCube Obser- 
vatory as well ()Abbasi et al.l 120 lit) . Although astro- 
physical arguments have been raised to expla i n the 
observed small scale features (Salvati ^Saccol 120081 : 
iDrurv fc Aharonianl[200llSalvati.2010.) . also m connec- 
tion to cosmic ray propagation t hrough the turbulent 
magnetic field (jMalkov et al.|[20To( ). the existence of such 
localized regions in the arrival direction distribution of 
TeV cosmic rays is a phenomenological indication that 
processes occurring within the mean free path are likely 
responsible for the observations. For instance, scatter- 
ing of TeV-PeV cosmic ray particles with the cascading 
turbulent LIMF within a few lO's pc can generate in- 
termediate and small scale perturbations over an under- 
lying global anisotropy (Giacinti & Sigl 2012). However 
the observed small scale anisotropy features at 1-10 TeV 

^ an Astronomical Unit (AU) is the mean distance between Earth 
and the Sun, corresponding to approximately 1.5x10^^ m. 



energy range, do not seem to have a stochastic nature, 
but they rather appear to be correlated with the global 
anisotropy structure. 

The change in topology of the cosmic ray anisotropy 
above about 100 TeV, is observed at the energy scale 
where the infiuence of the heliosphere is expected to be 
sub-dominant. Assuming that the anisotropy of TeV cos- 
mic rays approaching the heliosphere is the same as that 
above 100 TeV, this suggests that the arrival distribu- 
tion is globally re-shaped by heliospheric perturbations 
induced in the LIMF, and it might explain some appar- 
ently coincidental correlation between the anisotropy and 
the outer heliospheric features. In this paper, for the 
first time, combined full-sky observations of the global 
anisotropy of 1-10 TeV cosmic rays are used to provide 
some evidence that the change in topology at about 100 
TeV is due to the heliosphere, which is also responsible 
for the small scale angular structure. In Sec.[2]the helio- 
spheric magnetic field structure is illustrated. In Sec. [3] 
a more detailed description of the experimental observa- 
tions is provided and Sec. H] reports the discussion of the 
proposed model. 

2. THE HELIOSPHERE 

The solar system's motion through the LIC produces 
a complex interface due to the interaction between the 
supersonic solar wind and the interstellar flow, called 
the heliosphere. A termination shock, where the solar 
wind pressure equals that from the interstellar flow, is 
formed at approximately 100 AU from the Sun. The 
interface separating interplanetary and interstellar mag- 
netic fields, called heliopause, is at a distance of approx- 
imately 200 AU in the upstream direction (the helio- 
spheric nose) and it has a comet-like shape downstream 
due to the heliospheric plasma streamlines deviated by 
the interstellar flow (the heliospheric tail, or heliotail). 
The heliotail can be as wide as about 600 AU and as 
long as several thousands AU (jlzmodenov fc KallenbachI 
U006; Pogorclov et al. 2009a). Deformations on both the 
heliosphere and the LIMF drapi ng around it are pro- 
duced by the ir mutual interaction (jPogorelov et al.ll2004 
[2006ll2009bt ). The heliospheric magnetic field has been 
studied with detailed Magneto-Hydrodynamic (MHD) 
simulations, where the effects from the 26 day solar 
rota tion and the 11 year s olar cycle were considered 
(see IPogorelov et all ()2009a| )). Over the solar cycle the 
magnetic field polarity is reversed every 11 years, gen- 
erating unipolar regions dr agged along t he heliotail by 
the ^ 100 km/s solar wind (|Parkerill979l ). In particular 
these regions grow to their maximum latitudinal extent 
during solar minimum (about 200-300 AU in size) and 
reduce to zero at solar maximum, during which the he- 
liospheric plasma is dominated by the strongly mixed 
polarit y domains (about 0. 1-1 AU in size) from solar ro- 
tation (jNernev et al.l [19951 ) . Due to the tilt of the so- 
lar magnetic axis with respect to its rotation axis, the 
unipolar regions are thinner at lower latitudes. There- 
fore the tailward line of view is dominated by the finely 
alternating magnetic field, while along sightlines away 
from it the magnetic domains have larger sizc0. MHD 

^ the actual solar magnetic field structure and dynamics is signif- 
icantly more complex (Riley et al. 2011), thus affecting the detailed 
unipolar domains conformation, although the general properties 
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numerical simulations show that the sectored unipo- 
lar magnetic field regions can propagate along the tail 
for several solar cycles before they dissipate into the 
LISM. The periodic variations induced by solar cycles 
on the termination shock and the overall heliospheric 
plasma have profound effects on t he dynamic behavior 
of th e heUo sphere (Zank & Miille^ l2003l iWashimi et"an 
|2011LI2012| ). The corresponding magnitude of the LIMF 
can vary about 25% and th e Alfven velocity by about 
20% (jPogorelov et al.ll2009aD . 

There is observational evidence that the plasma in 
the heliosheath has Reynol ds number Re « 10^'* 
(see iLazarian fc Opheil (|2009l ) and references therein), 
meaning that the strength of non-linear convective pro- 
cesses at the largest scale is more important than the 
damping viscous processes in the dynamics of the flow. 
We expect a similarly high Reynolds number in the inner 
heliotail as well, and we know that this is true for the in- 
terstellar medium in general. In such conditions it is very 
unlikely that plasma flow stays laminar. In particular 
the wake downstream the interstellar flow develops tur- 
bulence at an injection scale of the order of the heliotail 
thickness (i.e. « 600 AU). In addition, the presence of 
neutral atoms in the partially ionized LISM is essential 
for the dynamics of the heliosphere and LIMF interac- 
tion. Charge- exchange processes between the interstellar 
inflowing neutral atoms and the outflowing solar wind 
protons can produce Rayleigh- Taylor type instabilities 
on the heliopause with amplitude of a few tens AU and 
over a time scale of a few hundreds years (jLiewer et al.l 
119961 : iZank et al1[l996lV In a model of plasma- neutral 
fluid coupled via collision and charge-exchange processes, 
it is found that such non-linear coupling leads to alter- 
nate growing and damping of Alfvenic, fast and slow 
turbulence modes, at L ~100's AU scale and with evo- 
lution time longer than inertial time L/Va, with Va 
the Alfven velo city (Zank 1999; Florinski et al. 2005; 
Borovikov erall[200& : .Zank et al.,,2009; .Shaikh fc Zanh 

The LIMF is deformed by the heliosphere over the 
largest scale of about 600 AU (i.e. the heliosphere 
transversal size), which is much smaller than the in- 
jection scale of the interstellar turbulence. The same 
non-linear processes that produce turbulence on the 
wake are responsible for its cascading down to smaller 
scales. Since turbulence evolves along hydrodynamic 
Kolmogorov cascades, the Alfven velocity 5Vi of helio- 
spheric magnetic perturbations at scale / can be ex- 
pressed as 



where Vh ~ 23-26 km/sec is the interstellar flow ve- 
locity (see iBzowski et al.l ()2012[ ) and references therein) 
and Lh ~ 600 AU the injection scale of the heliospheric 
induced turbulence. Since the Alfven velocity in the 
partially ionized plasma within the LIC that surrounds 
the h eliosphere is Va ~ 13-17 km/sec (Spangle r et al.l 
1201 1[ ). turbulence is super- Alfvenic at scales I a ^ 
Lh{VA/VhY w 80-250 AU. Magnetic perturbations at 
such scales are strong (i.e. SB/Bq w 1) and can- 
can still be described in terms of a tilted dipolar magnetic field of 
the Sun. 



not propagate into the upstream interstellar flow di- 
rection (i.e. they trail downstream). Super- Alfvenic 
turbulence isotropically cascades down to the scale I a, 
where it becomes sub-Alfvenic and the mean mag- 
netic field gets dynamically important. Cascading sub- 
Alfvenic turbulence is critically balanced, making per- 
turbations more elongated along the magnetic field 
at small spatial scal es, while keeping the ene rgy cas- 
cade rate constant (jGoldreich fc Sridhaii I1995D . This 
is true for incompressible turbulence, however it can 
be considered val id in the compressible case as well 
(see Lazarian fc Beresnvald (|2006[ )). which is what is 
expected to occur in realistic astrophysical plasmas. 
In this case turbulence can be composed by magne- 
tosonic modes, in addition to shear Alfven modes. On 
the other hand, turbulent dissipation in the cascad- 
ing process has the effect of suppressing magnetosonic 
waves and effectively makes turbulence nearly incom- 
pressible (Shaikh_fc Zank 200fi). As a consequence, com- 
pressible magnetosonic modes decay quickly leaving only 
the shear sub-Alfvenic modes to dominate the energy 
cascades. 

In the partially ionized very local interstellar medium, 
ion-neutral collisions is a source of damping mech- 
anism for Alfvenic turbulence, and the collisional 
mea n free path in the LISM is approximately 300 
AU (jSpangler et al.|[20Tll) . Therefore we can picture the 
outer heliospheric boundary as populated by cascading 
turbulent ripples in the LIMF trailing downstream (at 
the largest scales) and with exponentially decreasing am- 
plitude as a function of distance outward into the inter- 
stellar medium. 

3. OBSERVATIONS OF COSMIC RAY ANISOTROPY 

The observed global cosmic ray anisotropy from tens of 
GeV to tens of TeV energy range, is generally character- 
ized by a broad relative excess between right ascension 
a « 18 hr and a « 8 hr, and a relative deficit in the rest of 
the sky, as shown in Fig. [TJ The figure shows the map in 
equatorial coordinates, of cosmic ray relative intensity in 
arrival direction observed by Tibet AS 7 at median energy 
of 5 T eV in the northern hemisphere (jAmenomori et aLl 
l2Qlla[ ). and by IceCube at ab out 20 TeV niedian energy 
in the southern hemisphere (jAbbasi et al.l l2010t ) . The 
relative intensity measures the deviation of cosmic ray 
counts in a given declination band with respect to the 
average count in that band. Relative intensity distribu- 
tions are independently determined for each declination 
band (of order 1° — 5°) and are used to construct the 
map, thus providing information on the relative modula- 
tion of the arrival direction of cosmic rays along the right 
ascension. Although the right ascension position of the 
relative excess and deficit regions as a function of dec- 
lination depends on the cosmic rays energy, the global 
pattern of the anisotropy does not change significantly 
across the sky up to w 100 TeV. Moreover the energy 
response distribution for these observations are reported 
to be rather wide, therefore the different characteristic 
energies corresponding to the two halves of the sky are 
partially mitigated by a significant overlap in energy re- 
sponse. 

The relative excess region in the northern hemisphere 
shows a complex angular structure around the apparent 
direction of the interstellar downstream flow direction 
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Fig. 1. — Map in equatorial coordinates constructed by com- 
bining relative intensity distributions of cosmic ray counts 
independently normalized within declination bands of or- 
der 1° — 5°. The map shows the observation by Tibet 
AS7 at about 5 TeV i n the northern hemisphere (from 
lAmenomori et al] l|2011aD ) and by IceCube at a median en- 
ergy of 20 TeV in the southern hemisphere (from lAbbasi et al] 
((2010)). Note the different characteristic energy of the 
cosmic ray particles in the two hemispheres. Also shown 
the equator "at infinity" (see text) in black continu- 
o us line, co rresponding to the direction of LIMF inferred 
in lFunsten et al. (20 0^ (circul ar symbols), that in black dot- 
ted line from Schwa dron et al. (2009) ( square sym bols), and 
that in black dashed line from'Hccrikh uisen et al.l (2010) (tri- 
angle symbols). The star symbol indicates the downwind in- 
terstellar medium flow with coordinate (5 hr, -|-17°). Color 
scale indicates relative excess (red) and deflcit (blue) with re- 
spect to average intensity, and the map grid is 2 hr in right 
ascension and 15° in declination. 
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Fig. 2. — Schematic representation of the elongated helio- 
sphere with respect to the LIMF draping around it. The ar- 
rival distribution of cosmic rays with energy below 100 TeV 
have a relative excess in the portion of the sky south of the 
magnetic equator that includes the heliospheric tail in the 
downstream interstellar flow direction (represented in red). 
The deflcit is observed toward the upstream interstellar flow 
direction, north of the magnetic equator (represented in blue). 
The regions with the largest gradient in the cosmic ray inten- 
sity (in white between the excess and deficit) are approxi- 
mately located along the magnetic equator "at infinity" (see 
text). Also shown the celestial and ecliptic equators. 



(approximately the inferred direction of the heliospheric 
tail), with equatorial coordinates (a, d) « (5 hr, -1-17°) 
(indicated with a white star in the figure). This is the re- 
gion where a broad enhancement denominated tail-in ex- 
cess was obser ved for the firs t time with sub- TeV cosmic 
rays (Nagashi ma et allll998[ ). Although the global rel- 
ative excess observed in the southern hemisphere seems 
to have a different fine structure, it appears topologically 
connected to that observed in the north. This might 
be an indication that even though the energy distribu- 
tions of the two observations overlap, some small scale 
differences could still arise, perhaps also related to differ- 
ent mass composition sensitivity of the two experiments. 
The visible feature at equatorial coordinates {a, 6) ~ 
(21 hr, -1-32°) is reported to be originated by high energy 
gamma rays from the Cygnus region to which the Tibet 
AS7 is also sensitive. 

In general, the combined observations reported in 
Fig-Scan provide a valuable tool to study the global 1-10 
TeV cosmic ray anisotropy, although consistent observa- 
tions with comparable energy and mass sensitivities are 
desirable. It is interesting to note that such anisotropy 
structure is not compatible with a simple dipole, there- 
fore the cause cannot be simply from diffusive propaga- 
tion across the Galaxy. The figure shows the orientations 
of the LIMF as inferred by various interpretative models 
of the keV Energ etic Neutral Atoms (ENAs) observa- 
tions with IBEX (iFunsten et al.l 120091: iSchwadron et al.l 
[200I iHeerikhuisen et al.l 1201^ along with their cor- 



responding equator "at infinity'^ These orientations 
are consistent with that based on Ca II absorption 
lines and on obse rvations of H I lines in white dwarf 
stars (^rischlllQQGl) . with polarization data (jFrisch et alj 
l2012a||b(). and on observation of radio emissions from 
the inner helioshe ath and from observations of hydro- 
gen deflection tLalleinent et alJ 120051: iQpher et al.H2Q07l: 
iLallement et al.lboiOf) . If we assume that the various in- 
direct probes provide an accuracy of order 30° on the 
LIMF direction, we note that the magnetic equator is 
approximately located along a big portion of the rim, 
within about 8-10 hr right ascension, where the cosmic 
rays intensity exhibits a gradient within a narrow range 
in right ascension. Although this can be accidental, it is 
important to note that while it is not possible to decouple 
a relative modulation in declination (i.e. within right as- 
cension bands) from experimental coverage uncertainties 
in the experiment local zenith angle, the relative inten- 
sity in right ascension (i.e. within declination bands) is 
experimentally more robust and reliable. 

In order to relate the full-sky global anisotropy to the 
local environment, we schematically show in Fig. [5] the 
relative location of the excess and deficit regions of cos- 
mic rays with respect to the heliosphere and the LIMF 
that drapes around it. In the figure the termination 
shock (TS) is represented as a small circle at the cen- 

^ equator " at infinity" is the plane perpendicular to the magnetic 
field lines and passing by the observer's location, assuming the 
magnetic field is uniform. 
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ter and the heliopause is represented with a tail toward 
the left side. The heliotail axis (approximately the direc- 
tion of the interstellar wind upstream the heliosphere) , 
is in the plane of the figure, and it is coincident ally al- 
most on the ecliptic equatorial plane (which is perpen- 
dicular to the figure). In fact the ecliptic coordinates 
of the heliospheric nose are (A, /3) w (255°, 5°), i.e. 
just 5 degrees off the ecliptic plane. The LIMF is be- 
lieved to be approximately directed toward ecliptic co- 
ordi nates (A, /?) » (221 °, 39°) in the upstream direc- 
tion (|Funsten et al.|[2009[ ). This is about 40 degrees from 
the ecliptic plane and approximately 30 degrees above 
the plane of the figure, although the uncertainty from 
different determinations is at least of order 30°. For 
simplicity Fig. ^ schematically shows the component of 
LIMF in the plane of the figure with its draping around 
the heliosphere. The full-sky global anisotropy of the 1- 
10 TeV cosmic rays can be described as composed of a 
broad relative excess region across the hemisphere south 
of the magnetic equator (the red shaded sector in the 
figure, which contains the heliotail and the downstream 
direction of the LIMF), and a relative deficit region (the 
blue shaded sector) north of the magnetic equator. The 
cosmic ray intensity modulation in these regions cannot 
be described with a dipole, since the transition between 
the two regions lays approximately along the magnetic 
equator (the non-colored sectors in the figure) but with 
the steepest variations in directions with smaller angu- 
lar distance from the heliotail. At energies in excess of 
about 100 TeV the heliospheric infiuence must be sub- 
dominant, if not inexistent. Whatever the origin of the 
anisotropy above 100 TeV is, the arrival distribution of 
cosmic rays at lower energy is likely affected by the he- 
liosphere since the gyro-radius of « 8,000 AU is of the 
order of heliotail size. 

As mentioned above, the global TeV cosmic ray 
anisotropy appears to have a fine angular structure com- 
posed of localized features. Subtracting, from the obser- 
vation, a map where counts are averaged in each dec- 
lination band over a right ascension range of Aa, it is 
possible to highlight those regions where the cosmic ray 
intensity has variations on smaller scales than Aa. In or- 
der to assess whether the residual map significantly devi- 
ates from random fluctuations of the averaged map, the 
statistical significance of the observed deviations in rel- 
at ive intensi t y is c alculated using the method described 
in iLi fc Mai ( 1983D . In the determination of statistical 
significance the declination dependence of experimental 
counts is taken into account, thus providing the most ro- 
bust deter mination in d e clinat ion bands with the largest 
count fsee lAbbasi et al.l (|2011[ )). Fig. [3] shows the map in 
equatorial coordinates of the statistical significance with 
respect to the averaged map. The northern hemisphere 
shows the observation by Milagro at a mean energy of 
about 1 TeV, where the averaging was do ne in a 2 hr 
right ascension range (corresponding to 30°) (jAbdo et al.l 
l2008f ). The southern hemisphere shows the observation 
reported by IceCube at a median energy of about 20 TeV, 
with the averaging performed over 4 hr in right ascension 
(corresponding to 60°). To obtain the two maps the ref- 
erence level and residual data counts in each ^1° pixel in 
the sky were added within a radius of 10°, which there- 
fore r epresents the smoothing scale (see lAbbasi et al.l 
(|2011f )). Compared to Fig. [1] this map only shows the 
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Fig. 3. — Map in equatorial coordinates of the statistical 
significance of the observed cosmic ray counts with respect 
to that averaged over a right ascension interval Aa = 2-4 
hr. Non-zero values measure how significantly the relative 
excess and deficit regions (in red and blue, respectively) devi- 
ate from random fluctuations (see text). The map shows the 
observation b y Milagro at about 1 TeV in the northern hemi- 
sphere (from lAbdo et al.i (2008)) and that by the IceCube 
Observatory at a median energy of 2 TeV in the southern 
hemisphere (from 'Abt>asi. 'eFalT pOnl 'l'l. The symbols are as 
in Fig. [T] The map grid is 2 hr in right ascension and 30° in 
declination. 



small angular scale anisotropy fe atures. As in Fig. [1] 
the orientations of the L IMF fromlFunsten et al.l |2009.) ; 
iSchwadron et all (|2009( ) ; iHeerikhuisen et al.l (|20ia ) along 
with their corresponding equator "at infinity" are also 
shown. 

The most significant localized excess region observed 
in the northern he misphere (the so-called region A 
m lAbdo et all (jlOOl)) is in the apparent direction of the 
heliotail, or the downstream interstellar flow direction, 
with coordinates {a, S) w (5 hr, +17°) (the star in the 
figure). Another significant localized region of multi-TeV 
cosmic rays obser ved in the no rthern hemisphere (the so- 
called region B in lAbdo et al.l (|2008)) appears to be elon- 
gated within right ascension a « 8 - 10 hr. The excess 
at about (a, S) w (21 hr, -1-32°), also visible in Fig. [1] is 
generated by gamma rays from the Cygnus region, and 
it is n ot related to hadronic cosmic rays (see lAbdo et"al] 
(120081 )). The southern hemisphere shows similar small 
scale anisotropy features, although with lower statistical 
signiflcance due to the smaller number of events utilized 
to compile the map. The most signiflcant and broadest 
excess and deficit regions observed in the southern hemi- 
sphere, at {a, S) w (7 hr, -45°) and (a, 6) « (15 hr, 
-45°), respectively (see Fig. |3]) coincide with the center 
of the corresponding global anisotropy regions in Fig. |TJ 
This indicates that cosmic ray intensity variation at the 
relative deficit and excess is sharper than 60° scale. As 
already noted above, the flne structure observed in the 
north at right ascension « 4 - 6 hr does not seem to have 
a continuity toward the southern sky, which might be a 
residual effect of the different characteristic energy of the 
two observations. On the other hand the southern hemi- 
sphere elongated feature along right ascension a « 8 - 10 
hr, is apparently connected to the Milagro region B in the 
northern hemisphere. We note that such a localized fea- 
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ture does not appear to be particularly significant at this 
point. However the two elongated localized regions seem 
to be approximately located at the edge where the broad 
relative excess of cosmic rays steeply transitions into the 
deficit region (see Fig[T]). This indicates that the cosmic 
ray intensity variation across this edge is steeper than 
30°-60°. In addition a circle drawn over such elongated 
excess regions would approximately overlap with the rel- 
ative excess observed in the southern hemisphere at right 
ascension a sa 22 - 23 hr, also located by the region of 
largest intensity gradient, although it is not statistically 
significant and it might be a coincidence. In any case, the 
elongated localized fractional excess region across a big 
portion of the sky, is approximately parallel and close to 
sightlines perpendicular to the LIMF direction inferred 
from several indirect determinations (see Figs. [T] and |3]). 

It is important to note that the anisotropy observations 
at small angular scales might be affected by possible arti- 
facts arising from the experimental techniques used. On 
the other hand, based on the current results we find the 
correlation of the region of steepest cosmic ray intensity 
variations with the interstellar magnetic equator in the 
1-10 TeV energy range, compelling, and as a consequence 
we could relate some of the observed localized fractional 
excess regions, at least, to processes occurring within or 
in the outer boundary of the heliosphere. 

4. PROBING THE OUTER HELIOSPHERE AND THE LIMF 

Cosmic rays transport in the interstellar magnetic 
field is diffusive as long as they propagate through dis- 
tances larger than the scattering mean free path, be- 
yond which particles lose memory of their initial direc- 
tion, and can propagate across magnetic field lines due 
to their stochastic wandering at a scale of 10-100 pc. In 
the energy range of 1-100 TeV cosmic rays have a scat- 
tering mean free path in the inter stellar medium of or- 
der 1-10 pc (|Yan fc Lazarian||2008f) . therefore at smaller 
distances the diffusion regime is broken as particles ve- 
locities are not de-correlated. This means that parti- 
cles mostly stream along the LIMF within the interstel- 
lar mean free path. We assume that below about 100 
TeV cosmic ray arrival distribution is similar to that at 
higher energy. Whatever causes the anisotropy observed 
at several lOO's TeV, there is no reason to believe that at 
lower energy it has disappeared or significantly changed 
by purely interstellar propagation. Therefore, at energies 
where the heliospheric infiuence becomes important, as 
anisotropic cosmic rays approach the heliosphere, they 
begin to feel the turbulent perturbations SB induced by 
the wake, with exponentially increasing amplitude the 
smaller the distance. At the largest scales, turbulence 
is super- Alfvenic, therefore it is strong {SB/Bq w 1, 
with Bq ~ 3fiG the regular magnetic field) and isotropic 
(see Sec. [2]). Particles propagating through the turbulent 
magnetic field are affected by scattering processes that 
change their pitch angle, with a rate that depends on 
(SB/Bo)'^. Therefore, scattering rate increases fast the 
closer to the heliosphere, leading to a progressively more 
uniform pitch angle distribution. 

Cosmic rays predominantly scatter on magnetic field 
fluctuations that have a size scale L ^ Vg, with 

^ ^ ^ 200 /1mG\ 
ZeB ~ Z [iTeVj [ B J 



the particle gyro-radius, where p± is the momentum per- 
pendicular to the magnetic field direction, E the parti- 
cle energy, and B = Bq + SB the total magnetic field 
of the perturbation. Considering that SB/Bq « 1, this 
resonance condition is satisfied for proton energies of ap- 
proximately 10 TeV at the largest scales (i.e. « 300-600 
AU). This energy scale effectively determines the separa- 
tion between two regimes, as discussed in the following 
subsections. In resonant scattering with isotropic tur- 
bulence, particles pitch angle changes randomly by an 
amount of order SB/Bq in each collision 0. Therefore, 
when the average effective deviation in pitch angle af- 
ter N collisions is large (i.e. particles have scattered 
approximately 90° from the initial direction, or equiva- 
lently -s/iV {SB/Bq) « 1), particles have lost memory of 
their initial condition. The effective mean free path, de- 
fined as the distance at which pitch angle randomization 
is realized, is then 

X^NL^r,[l^)\ (3) 

Therefore, resonant scattering in strong turbulence at 
the largest scales generates a re-organization of the cos- 
mic ray arrival direction distribution. However, due to 
the super-Alfvenic turbulence distribution in the outer 
heliosphere at the largest scales, and the relative angle 
between the LIMF and the heliotail, the effective scat- 
tering rate depends on the arrival direction of the cosmic 
rays, and pitch angle scrambling is more important to- 
ward directions with the largest scattering. This may 
explain the complex angular structure of the observed 
TeV anisotropy, which could still hold some information 
on the direction and amplitude of the initial underlying 
anisotropic distribution. 

4.1. Scattering at about 10 TeV 

Protons at ss 10 TeV have resonant scattering with 
the isotropic super-Alfvenic turbulence at the largest 
scales. Particles streaming along the LIMF from the 
downstream direction (i.e. from the lower left side in 
Fig. m are expected to experience higher scattering rate 
while they get closer, due to the angular proximity to 
the heliotail. Because of the increasingly strong mag- 
netic field perturbations, resonant scattering scrambles 
high energy cosmic rays with the highest rate occurring 
at smaller impact distance from the heliotail. Cosmic 
rays that pass into the upstream region have very little 
chance to back-scatter, because the super-Alfvenic tur- 
bulent perturbations cannot propagate in that direction. 
In particular at the largest scales, turbulence trails down- 
stream at a mean angle of approximately 40° from the 
interstellar flow direction. For the same reason, energetic 
particles streaming along the LIMF from upstream, ex- 
perience very low scattering rate in the outer heliosphere 
until they reach the downstream region, where they can 
resonantly scatter and scramble their distribution, feed- 
ing and enhancing the downstream fiux. This is illus- 
trated in Fig. |4l where two energetic particles streaming 
along the LIMF from the nose side of the heliotail, and 
affected by the turbulent magnetic field, scatter at the 

which is the average inclination of the field lines from the mean 
field direction due to the magnetic irregularities. 
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Fig. 4. — Illustration of the turbulent ripples at large scale 
on the outer heliospheric boundary. High energy cosmic rays 
propagating from the upstream direction into the front of 
the heliosphere resonantly scatter with the trailing magnetic 
field perturbations on the flanks of the heliosphere and along 
the heliotail. The uniform pitch angle distribution makes 
it possible for particles to scatter back to Earth. The grey 
blurred area schematically represents the possible scattered 
particle trajectories toward Earth. 

flanks of the heliosphere and along the heliotail back to 
Earth. The uniform pitch angle scattering distribution 
scrambles particle directions and the possible trajecto- 
ries toward Earth are represented by the blurred grey ar- 
eas. If the turbulence surrounding the flanks of the outer 
heliosphere and the heliotail extends into the LISM for 
about 100 AU, resonant scattering can generate a flux of 
back-scattered particles that may be responsible for the 
re-distribution of the initial anisotropy, and in particular 
for the excess toward the heliotail. A quantitative esti- 
mate of back-scattered particles compared to the direct 
flux is done in Sec. 14.1.11 using a simplified model. This 
model shows that the overall effect of scattering on the 
turbulent outer heliosphere can be that of generating the 
observed asymmetry in arrival distribution. Because of 
the angle between the LIME and the heliotail and of the 
trailing of turbulence at the largest scales, the gradient in 
the cosmic ray intensity is approximately located toward 
sightlines that are perpendicular to the mean LIME. Such 
gradient is expected to be larger where the LIME hit the 
heliosphere at a steep angle (i.e. the upper left side in 
Eig. [2] and HJ corresponding to right ascension range of 
a ~ 8-10 hr in Eig.[T]), while the intensity has a smoother 
variation where the LIME tangentially drapes the helio- 
sphere (for instance at the lower right side in the figures, 
at a w 18-22 hr) because only smaller and weaker tur- 
bulent perturbations can effectively scatter cosmic rays. 

The line-of-sights normal to the mean LIME, can be 
several degrees away from the magnetic equator " at infin- 
ity" , because of the LIME deformation due to the helio- 
sphere. In fact, the locus of largest gradient in the cosmic 
ray arrival direction distribution, appears to be within 
about 30° from the magnetic equator inferred from in- 
direct estimations of the LIME direction. In summary, 
the relative angle between the LIME and the heliotail, 
although still uncertain, and the trailing turbulent mag- 
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Fig. 5. — Schematic model of cosmic ray scattering within a 
shell of turbulent plasma in the outer heliosphere. Energetic 
particles with guiding center streaming along the interstel- 
lar flow direction scatter inside the turbulent shell back to 
Earth. Due to resonant scattering and the length of the he- 
liotail, back-scattered particles can be of the same order of 
magnitude of the direct flux, or even slightly in excess. 

netic perturbations at the largest scale, are responsible 
for the non-dipolar nature of the global anisotropy. 

As mentioned in Sec. [U a steep intensity variation in 
the cosmic ray arrival direction distribution, may appear 
like localized fractional excess regions when a residual 
map is derived. The width of such localized regions is a 
function of the magnitude of the gradient, and it depends 
on the averaging procedure of the cosmic ray map. How- 
ever, in this model, it is the existence of such an intensity 
gradient that has a physical relation to the scattering 
processes of cosmic rays with the super-Alfvenic mag- 
netic field turbulence on the outer heliospheric boundary. 
The observation could have some relation with the rib- 
bon of keV ENAs observed by IBEX, attributed to some 
mechanism that enhances ENA flux in the sight of lines 
approximately perpendicular to the L IME as it drapes 
around the heliosphere. In particular, ISchwadron et aLl 
(2012) considered cosmic ray particles diffusion through 
the magnetic field draping around the heliosphere as 
the cause of the observed anisotropy, within a consis- 
tent model that addresses the IBEX ribbon. The more 
ecliptic southward location of the TeV cosmic ray obser- 
vation might indicate that its formation is farther away 
from the heliosphere, where scattering processes become 
important. 

4.1.1. Effect 0^ particle back-scattering 

Here, a simple model of the effect of cosmic ray back- 
scattering on the turbulent outer heliosphere along the 
tail is presented. Eor simplicity we assume that the guid- 
ing centers of the cosmic rays stream along the interstel- 
lar fiow (i.e. along the z-axis in Eig. [SI with its origin cor- 
responding to the location of Earth). In this model we as- 
sume that resonant scattering happens only within a shell 
around the heliotail of thickness (IRh (the grey shaded 
area in the figure), where the plasma is fully turbulent. 
The fiux of cosmic rays can be separated in an isotropic 
and an anisotropic contribution. Since the isotropic com- 
ponent remains as such, it is the anisotropic part of the 
flux that gets re-distributed by scattering processes. 

The anisotropic component of the cosmic rays, with the 
relative excess approximately toward the direction of the 
heliospheric nose (i.e. in the up-stream interstellar flow 
direction), is represented with flux density ncR- When 
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these particles pass through the turbulent heliospheric 
shell, they can have resonant scattering (depending on 
their energy and the turbulence scale). In particular, 
particles that reach an infinitesimal volume dV = dr dl dz 
of turbulent plasma past the Earth along the tail, are 
uniformly re-distributed and are back-scattered to Earth. 

The number of cosmic ray particles inside dV at a given 
time is 



dN = ncR dr dl dz. 



(4) 



If Pg is the scattering probability, the number of particles 
uniformly scattered from the infinitesimal volume dV is 



dN, = dN. 



(5) 



Therefore, the number of particles that reach Earth from 
dV is 



dNb=dNs^, 



(6) 



where Re ^ 4.3 x 10 ^ AU is the Earth's radius, and d — 
_|_ J.2 jg ^jj^g distance of dV from Earth. Combining 
Eqs. m Eland [S] we obtain 



dm = ncR Ps R% 



dr dl dz 
z^ -|_ 7.2 



(7) 



Since we want to calculate the amount of cosmic rays 
that scatter back to Earth, we integrate Eq. [7] over the 
volume of the turbulence shell of radius and thickness 
dRii and along positive values of z (i.e. past the Earth, 
see Fig. [5]). So the number of cosmic rays back-scattered 
to Earth is 



m = ncRPsRl 



dr 



dl 



dz 



_|_ 



= ncR Ps TT^ i?| dRH , 



(8) 



where the heliotail is long enough (about 10^-10'* AU) 
that we can approximate the integration on z to infinity. 
Eq. [8] implicitly assumes that the contribution of cosmic 
rays that back-scatter from the far regions of the heliotail 
is taken into account. From the integral in z in Eq. [5]we 
see that about 50% of particles back-scattered to Earth 
comes from within a distance of about 300 AU along 
the heliotail. Cosmic rays reaching that distance and 
propagating back to Earth take about r ~ 6 x 10^ sec, 
which is considered as the back-scattering time scale. In 
this time the number of cosmic rays directly reaching 
Earth (i.e. from the negative values of z) is 



Nd = ticrAt^ReCt. 



(9) 



The asymmetry derived from the back-scattered par- 
ticles with respect to the direct flux can be evaluated 
with 



Nb-Nd ^ Nb/Nd - 1 
Nb + Nd Nb/Nd + V 



(10) 



TVb ^ Stt dRn^ 
Nd~ A ' CT 



(11) 



Based on this model, the back-scattered cosmic rays 
exceeds the direct flux (i.e. S >0, meaning that the ini- 
tial anisotropy is completely re-directed approximately 
toward the heliotail) if Ps }Z 100 /dRn- Therefore if we 
assume that dRn ^ 100 AU, the scattering probability 
must be very close to 100%. Although it is reasonable 
to assume that at the resonance Ps ~ 100%, experimen- 
tal observations are affected by energy resolutions that 
are not negligible, therefore they cover ranges where the 
effects of resonant scattering are reduced, thus decreas- 
ing the effective scattering probability. Therefore, based 
on this model, in order to have enough back-scattered 
particle to re-organize the cosmic ray arrival distribution 
at a given energy, the turbulent shell should be wider 
than 100 AU. However, the energy resolution of the ob- 
servations relaxes such condition by allowing particles at 
different energies to contribute to the back-scattered flux. 

We believe that the simple assumptions made in this 
model, do not affect the general evaluation of the ef- 
fects of cosmic ray resonant scattering with the turbu- 
lent outer heliosphere along the tail. On the other hand, 
the complexity of the heliospheric structure and of the 
turbulence on its outer boundary, and the fact that cos- 
mic rays mainly stream along the interstellar magnetic 
field, which is at an angle with respect to the heliotail, 
requires full numerical simulations in order to provide a 
solid prediction of the effects that scattering processes 
have on TeV cosmic rays. 

4.2. Scattering below 10 TeV 

Cosmic rays in the 1-10 TeV energy range can have 
resonant scattering with the cascading isotropic turbu- 
lence at scales down to order of 100 AU. The decrease 
in power of the smaller turbulence scales is compensated 
by the higher cosmic ray density at lower energy, so that 
effectively resonant scattering is still significant, and the 
anisotropy has the same topology and similar amplitude. 

Sub- TeV cosmic rays can have resonant scattering only 
at scales where turbulence is sub-Alfvenic. Since at such 
small scales cascading turbulence does not have time to 
develop a significant anisotropy, resonant scattering at 
scales of order 10 AU is a possible important contribu- 
tion. However damping processes might be large enough 
to limit turbulence inertial range at this scale. On the 
other hand an important contribution comes from par- 
ticle propagation through larger turbulence scales. In 
this case particles trajectory follow their guiding cen- 
ter along magnetic field lines and the pitch angle varies 
gradually with the changes in magnetic field due to the 
conservation of the first adiabatic invariant p\/B. An 
increase of B over times longer than particle's gyration 
time, produces an increase in the transversal momen- 
tum, leading to an effective change in pitch angle. In 
this case pitch angle scattering is not completely uni- 
form, leading to a similar anisotropy topology, but with 
reduced amplitude, as confi rmed by observations (see for 
instance lAmenomori et ahl 1 



where 



The mechanism where relatively low energy cosmic 
rays are affected by larger turbulence scales, where non- 
resonant scattering can be important in scrambling cos- 
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mic ra y arrival directions ( see also [S chlickciscr & Miller 
(I1998D ) , was considered b y lYan fc Lazarian (.2004 .2008, 
1201 ID : iBeresnvak et alj ()201lD . In particular, such 
non-resonant processes with large turbulence scale pro- 
duce gyro-resonance instabilities that effectively transfer 
turbulence energy into additional small scale Alfvenic 
waves, th us increasing low energy co s mic ray scatter- 
ing r ate (jLazarian fc BeresnvakI l2006t lYan fc LazarianI 
1201 It) . This mechanism might not play the same role 
in the case of small inertial range, although this would 
require a more dedicated study with numerical simula- 
tions. In any case, the fact that anisotropy observations 
below 1 TeV provide evidence of a tail-in excess, supports 
the idea that scattering processes are likely important in 
the sub- TeV energy range as well. 

4.3. Scattering above 10 TeV 

At energies above about 10 TeV, proton gyro-radius 
is greater than the largest turbulence scale. In this case 
particles move in orbits determined by the mean mag- 
netic field, and the resonant scattering with magnetic 
fluctuations starts to be sub-dominant. Lower scram- 
bling of cosmic rays pitch angle means a decreasing 
anisotropy, as observed. In particular, the intensity gra- 
dient by the LIMF equator becomes progressively more 
gradual and, eventually, the small scale features broader 
and indistinguishable from the global anisotropy dis- 
tribution. Non-resonant scattering with turbulence at 
scales much smaller than the gyro-radius can have im- 
portant effects as well, independently of the degree of 
turbulence anisotropy, although this case has not been 
studied in detail yet. 

For energies higher than about 100 TeV, even the pro- 
ton gyro-radius is larger than the heliotail and cosmic 
ray anisotropy is dominated by interstellar effects within 
distance scales of their mean free path in the interstellar 
medium. 

Heavier cosmic rays can contribute to resonant scat- 
tering processes at higher energies, possibly producing 
significant effects from heliospheric perturbations well 
above 100 TeV, depending on mass composition. For 
instance cosmic rays are expected to be dominated by 
helium nuclei above about 10 TeV ([Ahn et al.l 120101 : 
lAdriani et al.ll2011l : IPTSciascio et al1l2012t ). and to have 
significant contribution from he avier particles at w PeV 
range (i.e. the knee region) (jHorandell I2003D . The 
evolution to heavier mass composition, inevitably has 
the effect to smear the energy transition between he- 
liospheric ajid_2iite;stellarajuB hund reds 
TeV (see lAglietta et al.l (|200i ): lAbbasi et al.l (|20T1 ). 

Observations of cosmic ray anisotropy with experi- 
ments that have a wide energy response, result in averag- 
ing scattering contributions at different scales, and there- 
fore in smearing the resonant scattering effect. The re- 
sulting sky maps inevitably show anisotropies originated 
at different energies and make it difficult to disentan- 
gle the individual effects. However, current observations 
make it still possible to provide good order of magnitude 
comparisons with models. 

4.4. Scattering in the inner heliotail 

Cosmic rays in the GeV to TeV energy range, 
propagating through the inner heliotail, are also af- 
fected by the complex magnetic field structure, with 



^200 AU unipolar domains generated by solar cy- 
cles, and strongly mixed magnetic field regions in be- 
tween. Stochastic magnetic reconnection processes may 
cause re-acceleration as long as particles can be trapped 
inside the reconnecting domains, and scattering pro- 
cesses are not sufficiently strong to cause significant 
energy re-distribut i on (se e Lazaria n fc Desia ti (2010); 
iDesiati fc LazarianI ()2012j )). Such scattering processes 
are weaker for those energetic particles propagating at 
low solar latitudes, where the magnetic field is strongly 
mixed at small scale (see Sec. [2]). Away from this direc- 
tions scattering with the larger unipolar domains scram- 
ble TeV particles in a broader distribution, providing an 
environment where acceleration by magnetic reconnec- 
tion is less effective. This may be at the origin of the lo- 
calized fractional excess region observed in the interstel- 
lar flow do wnstream direction and of the claimed sp ectral 
anomaly (Abd o et al1 l2008': 'Pi S ciascio et al.l l2012t 

Observations of a correlation of low energy cosmic rays 
with the solar cycles, have been documented and un- 
derstood as due to the influence of the variable mag- 
netic field intensity and turbulence from solar activity. 
An anisotropic effect correlated with solar cycles was 
observed with 600 GeV cosmic rays as well, and inter- 
preted in terms of acceleration due to the solar wind 
induced electric field expected when c osmic rays cross 
the solar magnetic field neutral sheets (iMunakata et al.l 
l2010f ). Recently, a ~ 4xl0~^ relative modulation of 
TeV muon intensity (produced by ~ 10 TeV cosmic 
rays) correlated with the solar cycle, superimposed to 
the well known seasonal variation due to the upper at- 
mosphere's temperature, was reported based on over 20 
years o f muon measure ments at the Gran Sass o Labo- 
ratory ([Fernand ez-Mart inez fc Mahbubanill2012| ) . If this 
observations will be confirmed by other long term mea- 
surements of multi-TeV cosmic rays, it might provide a 
confirmation that at such energies particles from the di- 
rection of the heliotail may actually be influenced by the 
solar magnetic field polarity domains within a distance 
of about 300 AU. 

5. CONCLUSIONS 

TeV cosmic ray anisotropy can be the result of ac- 
cidental galactic sources that happened to accelerate 
particles at relatively close distance in relatively recent 
times. This intriguing scenario that relates the observed 
anisotropy to the origin of galactic cosmic rays can be 
used to provide constraints on particle acceleration at the 
source and propagation across the galactic disk, if diffu- 
sive regime is assumed. These models should be able to 
consistently describe the diffuse 7 ray emissions from the 
interactions of cosmic rays with the interstellar medium, 
molecular gas, dust and radiation and magnetic fields, as 
well. On the other hand the uncertainties on the source 
distribution, and the simple assumption of single power 
law scattering frequency, and on their dependency on the 
type of galactic medium, make such predictions suscep- 
tible of large uncertainties. The non-dipolar structure 
of the global cosmic ray anisotropy and its fine angular 
structure challenges theories relating this observations 
with large scale diffusive models. Even though it is likely 
that one or more nearby sources imprint an anisotropy, 
non diffusive cosmic ray propagation at distances smaller 
than their mean free path can easily re-distribute their 
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arrival direction. Within such a propagation volume, 
stochastic fluctuations in the interstellar medium mag- 
netic field can generate random fluctuations in the cosmic 
ray anisotropy that appear as localized excess regions. 

While awaiting for more refined experimental obser- 
vations of cosmic ray anisotropy as a function of energy 
and angular structure from various experiments and with 
higher event statistics, it seems natural to assume that 
the elongated heliosphere must affect in some way cosmic 
rays with gyro-radius comparable to its size. Cascading 
turbulence from the wake, which trails downstream the 
interstellar flow and outward into the LISM, resonantly 
scatter TeV cosmic rays streaming along the LIMF. The 
angle between the heliotail and the LIMF and the asym- 
metry in the turbulence from the up-stream and down- 
stream direction produces the non-dipolar anisotropy 
observed below about 100 TeV. At higher energy the 
anisotropy has a different structure and is most probably 
affected by propagation processes at larger scales, or even 
have an imprint from one or more galactic sources. The 
structure of TeV cosmic rays is related to the interplay 
between the LIMF and the heliosphere, and therefore it 
can be used to infer information on the interstellar mag- 
netic fleld direction just outside the heliospheric bound- 
ary (i.e. within several hundreds AU). The small scale 
features that arise where the gradient of cosmic ray in- 
tensity is sufficiently large within 30°-60° angular range, 
highlight a possible connection to scattering with the tur- 
bulent magnetic field. 

In order to validate such model, dedicated numerical 
MHD simulation models of the interaction between the 
interstellar and heliospheric magnetic fleld, and of the 
their turbulence, need to be developed. Current MHD 
simulations of the heliosphe re have reached rather high 
complexity (see for instance iPogorelov et al.l ()2009al lb[) ) . 
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Although in general it is not trivial to combine ob- 
servations from various experiments to complete the sky 
coverage, because of the different energy response and 
analysis techniques, it is possible now to have a first view 
into the global anisotropy and its small angular scale 
structure in the 1-10 TeV range. Dedicated observations 
with comparable energies from different parts of the sky 
can provide more solid experimental evidence to test the 
models. In the meantime we consider the model pre- 
sented here as a simple explanation of the observations, 
that will need verification as progresses in heliospheric 
modeling and more experimental results are available. 
The simple quantitative model presented in Sec. l4.1.T] is 
just able to provide a glance into the proposed mecha- 
nism. Due to the complexity of the heliospheric structure 
and the turbulence that it entails, only numerical simu- 
lations can provide a solid prediction of the effects that 
scattering processes have on TeV cosmic rays. 
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